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ABSTRACT 

We present cuspy, non-axisymmetric, scale-free mass models of discs, whose grav- 
itational potentials are of Stackel form in parabolic coordinates. A black hole may 
be added at the centre, without in any way affecting the Stackel form; the dynamics 
in these potentials is, of course, fully integrable. The surface density, Edisc oc 1/r 7 , 
where < 7 < 1 corresponds to steep cusps for which the central force diverges. Thus 
cusps, black holes, and non-axisymmetry are not a sure recipe for chaos, as is gener- 
ally assumed. A new family of orbits, lens orbits, emerges to replace the box orbits 
of models of elliptical galaxies that have constant-density cores. Loop orbits are con- 
spicuous by their absence. Both lenses and boxlets (the other family of orbits), can be 
elongated in the direction of the density distribution, a property that is favourable for 
the construction of non-axisymmetric, self-consistent equilibrium models of elliptical 
galaxies. 
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1 INTRODUCTION 

More than two decades ago, it became apparent that el- 
liptical galaxies are not rotationally supported, oblate ob- 
jects (Bertola & Cappacioli 1975, Illingworth 1977). Sub- 
sequent dynamical investigations (c.f. Binney 1976, 1978, 
Schwarzschild 1979, 1982) led to the picture of ellipticals 
as slowly/non-rotating, triaxial objects with anisotropic ve- 
locity distribution, whose self-consistent gravitational po- 
tentials supported mostly regular orbits. The discovery of 
fully integrable, triaxial mass models by Kuzmin (1973) and 
de Zeeuw & Lynden-Bell (1985) provided analytic power, 
and opened up a new era of dynamical investigations. The 
key insight was the employment of "Stackel potentials" in 
ellipsoidal coordinates; in addition to the energy, stellar or- 
bits in Stackel potentials respect two extra isolating inte- 
grals. A range of triaxial, self-consistent equilibrium con- 
figurations could now be constructed in an efficient manner 
(Statler 1987); a property common to all these models is a 
constant-density core in which stars execute near-harmonic 
oscillations ( "box orbits" ; see de Zeeuw 1985 for a thorough 
account of orbits). 

Recent observations (c.f. M0ller, Stiavelli, & Zeilinger 



1995, Ferrarese et al 1994, Lauer et al 1995) of elliptical 
galaxies, however, do not support the notion of a constant- 
density core; indeed, the density of stars appears to rise 
toward the centre in a power-law cusp. It is also widely 
believed that these galaxies harbour dark objects, possibly 
supermassive black holes, at their centres (c.f. Kormendy & 
Richstone 1995). Gerhard & Binney (1985) showed that cen- 
tral black holes, or steep density cusps will destroy box orbits 
in triaxial configurations, and argued that the boxes will be 
replaced by chaotic orbits. They also noted that chaotic or- 
bits fill a region of space that is too round for strongly non- 
axisymmetric, self-consistent equilibria to be constructed; 
thus it has been suggested that cuspy galaxies with black 
holes must necessarily be axisymmetric, at least near the 
centre (c.f. Merritt 1996 for an account based on recent 
numerical explorations). Although the link between cusps, 
non-axisymmetry, and chaos is supported by numerical re- 
sults in model potentials, we question its generality. We con- 
struct separable mass models of cuspy, non-axisymmetric, 
scale-free discs, with central black holes in this paper. The 
potentials of these discs promote a new family of orbits, lens 
orbits, to replace the box orbits of galaxies with constant- 
density cores. We hope that our results will add fuel to the 
ongoing debate on the intrinsic shapes of elliptical galaxies. 
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2 POTENTIAL-SURFACE DENSITY PAIRS 

We begin with an explicit formula for the potential (in the 
plane) of our discs, in the usual (r, 9) polar coordinates: 



$ = Kr 



where K > is a constant, and 

F{9) = {{! + sin0) (2 ~ 7) + (1 - sm<9) (2 ~ 7) } 



(1) 



(2) 



Knowledge of the potential in the plane of the disc suffices 
to determine the surface density though the integral (c.f. 
Binney & Tremaine 1987), 



E(r',9') 



(V 2 $) rdrd.9 
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(3) 



It is evident, from dimensional considerations, that E(r, 9) = 
5 , (6')/r 7 + M5(r) is the form of the surface density arising 
from a $ given by equation (|l|). For < 7 < 1, straightfor- 
ward manipulations of equation (|^), using Legendre func- 
tions P—y, demonstrate thalj^| 



s(e') = 



K 



de{(i-j) 2 F(e)+ 



AtyG sin(7T7) 

f"(q) }P-y[-cos(e -e')] ; 



(4) 
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de{sm 6>P" 7 (- cos #) + 



8ttG sin(7T7) 

(i-j)(2-j)p^(-co S 9)}F(e + e'), 

is finite and positive. Finiteness is guaranteed by the singu- 
larity of P Y ((j,) at (j, = — 1 being integrable, and positivity 

simply because F(9), P_ 7 (/x) and P" 7 (/z) are all positive 
quantities. If we neglect the contribution from the black hole, 
< 7 < 1 for all our models implies that the force on a test 
particle diverges, whereas its circular speed goes to zero for 
small r . Isocontours of surface density, and potential for 7 
equal to 0.1 and 0.5 are displayed in Figure 1. The contours 
(Figures la and lc) are guitar-shaped, with the depression 
near the ordinate increasing with 7; of course, the potential 
isocontours are much rounder than the density isocontours. 
The location of the black hole — the origin — is indicated by 
the solid dot. While drawing the potential isocontours (Fig- 
ures lb and Id), we have excluded the contribution of the 
black hole, to better display the potentials of the discs them- 
selves. The major axes of the potential isocontours are ex- 
actly twice as long as the minor axes, a property that is 
readily verified using equation (jij). 

A remarkable property of this family of potentials is 
that all stellar orbits are non chaotic. To make this explicit, 
we will rewrite it in Stackel form in parabolic coordinates. 



t 



The proof of the second of these relations requires integration 



by parts, and the elimination of PL-y by using the standard dif- 
ferential equation satisfied by P_ 7 (c.f. Morse & Feshbach 1953). 







Figure 1. Isocontours of Surface Density and Potential: Figures 
(a) and (c) show isocontours of the surface density for two differ- 
ent values of 7. The corresponding potentials are displayed in the 
panels on the right in (b) and (d). Successive isocontours in (a) 
have density ratios of 1.1, whereas the ratio is 1.5 for the other 
three figures. The location of the central black hole is shown as 
a solid dot in all four figures, although the contribution to the 
potentials is not included in (b) and (d). 



3 DYNAMICS IN PARABOLIC 
COORDINATES 

Let f = ^(sin^ + 1) and 77 = r(sin# — 1) be parabolic co- 
ordinates in the plane. Lines of constant £ and rj intersect 
the j/-axis at £/2 and n/2 respectively, and the x-axis at ±£ 
and ±17 respectively. Denoting the canonical momenta by 
and p v , the Hamiltonian, 



H(Pi,p v ,£,v) = 



2$ 



2 , 



(5) 



governs dynamics in the (£, 77) plane. Since <E> is independent 
of time, H — E is a conserved quantity. The potential of 
equation (Jlf) may be written in the Stackel form, 



To verify this, we choose 



(6) 



F+(0 = 2Ki 



(2-7) 



GM 



(2-7) 



+ GM. 



(7) 



and substitute for (£, n) in terms of (r, 9) . It is well-known 
that, for such potentials, the Hamilton-Jacobi equation sep- 
arates, yielding an extra conserved quantity (c.f. Landau & 
Lifshitz 1976), 



2frl-EZ + T + (0 
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Erj + T-(rj), 



(8) 



The isolating integrals, E and I may be used to classify 
orbits. For the generic case of a black hole, plus a cuspy den- 
sity distribution, E can assume both positive and negative 
values. For fixed E, the "second" integral, I, will determine 
the excursions in £ and n. The requirement that p\ and p\ 
be non negative, forces 



i>g(0, -i>g(\v\), 

where 

g(s) = 2Ks (2 ~"' ) -Es- GM : 



s > 0. 



(9) 



(10) 



For fixed E, the range of I is determined by the minimum 
of the function g(s), which is necessarily non-positive. If 
we denote this minimum value by —I lri (E), we obtain the 
condition — I m < / < 7 m . Since s > 0, 



Im(E) 



GM 



+ GM if E > 
if E < . 



(11) 



where C = (1 - j)[2K(2 - j)] 1/ ^- 1) /(2 - 7) is a positive 
constant. Let us consider two cases (the illustrative figures 
for which, given in Figures 2a-2d, are drawn for 7 = 0.5): 

(i) E < : We find that g(s) is a monotonically increasing 
function of s, attaining a minimum value of —GM for s = 
(see Figure 2a). This implies that —GM < I < GM. For 
a fixed (say, positive) value of /, the motion is bounded by 
the coordinate curves < £ < £ m , and < \n\ < rjm, where 
(,m(E,I) > f] m (E,I) > are the two roots of g(s) = ±1; 
the intersections of the dashed lines with the solid curve, in 
Figure 2a, gives the location of £ m and ?7 m . The orbits fill 
a lenticular region, bounded by the parabolas £ = £ m and 
Tj = — r\ m \ we call these lens orbits. As the shaded region 
of Figure 2b indicates, the lens orbits can visit the origin; 
they are to cusps with black holes, what box orbits are to 
analytic, triaxial cores. When |7| = GM, its maximum value, 
the shaded region collapses to an interval of the y-axis, < 
V < Cm/2. For negative values of I, the roles of £ and 77 are 
interchanged. 

(ii) E > : As seen in Figure 2c, g(s) is no longer a 
monotonic function of s; a minimum, — I m < 0, is attained 
at s = [E/2K(2 - 7)] 1/(1 " 7) . For |/j < GM, we again 
obtain lens orbits. But a new type of orbits appears for 
GM < \I\ < J m ; the shaded region of Figure 2d shows 
the region filled by one such orbit for a typical value of I, 
which is taken to be positive. These are centrophobic box- 
like orbits; we call these boxlets, borrowing the term from 
Miralda-Escude & Schwarzschild (1989). When I takes its 
maximum value equal to J m , the orbit collapses to one sec- 
tion of a constant -r) curve lying between two halves of the 
£ = Cm curve. This is a 2 : 1 resonant banana orbit, which 
now emerges naturally as a section of a parabola! 

It is worth noting that loop orbits are completely ab- 
sent. In the absence of the black hole {M = 0), the lens 
orbits disappear, leaving the boxlets as the only family of 
orbits for cusps without black holes. 
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Figure 2. Orbital Structure: In all figures 7 = 0.5, and / has been 
asumed to be positive. The solid vertical and horizontal lines in 
(a) and (c) correspond to the axes s = and g(s) = respec- 
tively. Similar lines in (b) and (d) are the axes x = and y = 
respectively, (a) E < 0, The solid curve is g(s), which intersects 
the dashed lines g = ±7 at s = 5m and s = n m respectively, (b) 
Lens orbits are restricted to the shaded region of the x — y plane, 
which is bounded by the curves £ = £ m and n = — rj m . (c) When 
E > 0, the minimum of g(s) is shifted away from s = 0. For 
GM < I < Ira, the dashed line at —I intersects the g(s) curve at 
two non-zero values of n. This property forces the boxlets to be 
centrophobic, as shown in (d). 



4 CONCLUSIONS 

The potential-density pairs presented in this paper are the 
only examples known to us, of non-axisymmetric distribu- 
tions of matter that have density cusps giving rise to sep- 
arable motion. The potentials are of Stackel form, and owe 
their existence to the geometric properties of the parabolic 
coordinate curves.^] As a bonus we can also include a central 
black hole of arbitrary mass, while maintaining the separable 
nature. The discs are significantly non-axisymmetric, with 
the major axes being more than twice as long as the minor 
axes. The cusps can be steep, with logarithmic slopes of the 
surface density lying between —1 and . Different discs may 
be superposed to obtain mass models that are not scale-free, 
yet whose potentials are of Stackel form. It is a surprising 
fact that the presence of a black hole stabilises a family of 
orbits, the lenses, which can approach arbitrarily close to the 
centrelpjj The lenses do not exist in the absence of the black 
hole — the centrophobic boxlets are the only stable family of 
orbits for non-axisymmetric cusps, without a central black 



3 Stackel potentials in elliptic coordinates do not appear to de- 
scribe density cusps of interest to dynamics at the centres of el- 
liptical galaxies! 

" Of course, a star that approaches too close to the black hole 
will be tidally shreded. 
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hole. From the geometry of the constant £ and r\ curves, 
it is evident that the lenses and boxlets can be elongated 
in the same sense as the mass distribution, a property that 
is favourable to the construction of non axisymmetric, self- 
consistent equilibria. Although we have limited discussion to 
discs that are symmetric with respect to reflections about 
both x and y axes, it is a simple matter to break the symme- 
try with respect to the x axis; weighting the functions J-+{^) 
and T- (77) differently in equations ((?]) will produce a family 
of lopsided, cuspy discs. In fact, choosing functional forms 
for and J--(rf) different from the ones given in equa- 

tions (pi) presents us with a variety of separable potentials 
(which, however, may not always be derivable from posi- 
tive densities). The separable potentials we have presented 
might also be a suitable point of departure for perturba- 
tive analyses of dynamics in more realistic models of non- 
axisymmetric, lopsided discs. The scale-free nature simpli- 
fies the problem, but we clearly need to construct more re- 
alistic, finite mass configurations. Generalization to triaxial 
cusps, and construction of self-consistent equilibria remain 
as some of the outstanding questions. The restoration of in- 
tegrability, and the existence of favourably oriented orbits in 
the idealized problem we have considered suggests that it is 
perhaps premature to conclude that elliptical galaxies with 
central density cusps, and black holes cannot be triaxial. 
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